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Abstract. We describe light scattering from a graphene sheet having a modulated 
optical conductivity. We show that such modulation enables the excitation of 
surface plasmon-polaritons by an electromagnetic wave impinging at normal incidence. 
The resulting surface plasmon-polaritons are responsible for a substantial increase 
of electromagnetic radiation absorption by the graphene sheet. The origin of the 
modulation can be due cither to a periodic strain field or to adatoms (or absorbed 
molecules) with a modulated adsorption profile. 
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1. Introduction 

Since the days of A. Sommerfeld, back to 1899, tliat plasmonics effects in materials and 
gratings are investigated both theoretically and experimentally. Modern plasmonics 
gained a renewed interest with the discovery super transmittance through a periodic 
array of holes with a size smaller than the diffraction limit [1]. The effect was explained 
invoking surface plasmons (SPs) [21 [Sj IH [5] . Beyond fundamental physics, plasmonics 
encompasses wide range of applications, such as spectroscopy and sensing [HI [3 E], 
photovoltaics [9], optical tweezers [lOllTT], nano-photonics [HllS], transformation optics 
[13], etc. 

A series of recent papers [TH|15l|16] triggered a burst of interest on plasmonic effects 
in graphene. In particular, it has been shown that graphene has a strong plasmonic 
response in the THz frequency range at room temperature [15]. THz photonics is 
emerging as an active field of research [17] and graphene may play a key role in THz 
metamaterials in the near future. Ju et al. [15J has shown that electromagnetic radiation 
impinging on a grid of graphene micro-ribbons can excite SPs on graphene leading to 
prominent absorption peaks, whose position can be tuned by doping. In general, SPs 
cannot be excited by directly shining light in a homogeneous system due to kinematic 
reasons: the momentum of a SP is much larger than that of the incoming light having 
the same frequency. Therefore some type of mechanism is necessary to promote the 
excitation of surface plasmons. The most common mechanisms for SP excitation are: 
attenuated total reflection (ATR) [TB] , scattering from a topological defect at the surface 
of the metal [3] , and Bragg scattering using diffraction gratings (or producing a periodic 
corrugation) on the surface of the conductor \V5]. The method of Ju et al. is similar 
(but not exactly) to patterning a metallic grating [16] on top of graphene. A theoretical 
account of the experiment by Ju et al. [15] was given in a recent work [20] . 

The reason why a periodic corrugation allows the excitation of SPs can be 
understood by an analogy with the theory of electrons in a periodic potential. The 
corrugation plays the role of the periodic potential. Therefore, the SP momentum 
is conserved up to a reciprocal lattice vector, that is to say, the periodic corrugation 
provides the missing momentum needed to excite the SP. Another way of seeing the 
effect is to note that the grating gives rise to a SP band structure in the first Brillouin 
zone. Then, the folding to plasmonic bands makes it possible to the impinging light to 
excite a SP associated with the upper bands in the Brilloun zone. The excitation of 
SP modes in the first band is not possible in the grating configuration unless the ATR 
technique is used. 

The question we want to answer in this article is the following: in what 
circumstances can electromagnetic radiation (ER) impinging at the interface of a 
dielectric and a metal excite SPs when the surface of the conductor has no corrugation? 
Taking the example of graphene, we show that a periodic modulation of the conductivity 
suffices for that end. In situations where both modulation of the conductivity and 
corrugation are present, understanding the former situation alone is a needed research 



Light scattering by a medium with a spatially modulated optical conductivity: the case of grapheneS 



step. 



2. Modulated conductivity 



Several procedures can be used to induce a patterned conductivity on a graphene sheet. 
One way is using split gates [21]. This leads to a modulation in the electronic density 
which induces a modulation of the conductivity. One can also assume that a CVD grown 
graphene sheet is transferred onto a patterned substrate, as illustrated in Fig. [l] The 
transfer process together with the patterning act in a way to produce inhomogeneous 
strain in the graphene sheet. Qualitatively, we expect the strain to be higher in regions 
labeled by 1 in Fig. [T] than in those labeled by 2. In a typical experiment we have 
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Figure 1. Periodic patterned substrate (top) and adsorbed atoms or molecules 
(bottom). For the patterned substrate we assume h <^ D and, of course, the graphene 
sheet will not follow the profile of the gate exactly. The picture is meant to represent 
a possible away of inducing a periodic strain field in graphene; the experimental 
realization may be quite different from this representation. In the bottom figure we 
consider a different realization: the graphene sheet is locally doped by adsorbed atoms 
or molecules which create an inhomogeneous charge density profile and therefore a 
modulated conductivity. 

h/D 10^^ ^ 1. In this regime we can safely neglect the effect of the groves in what 
concerns electromagnetic radiation scattering. Due to strain, the optical conductivity 
of graphene will be spatially modulated with the period of the patterned substrate. 
Naturally, the form of the modulated conductivity depends on the strain field. It is a 
rather complex problem to determine the exact form of the spatial conductivity, mainly 
because the strain field is not known exactly. In any case, we can say that the strain 
field should have the same period of the patterned substrate. 

Another possibility for producing a modulated conductivity is by molecular 
adsorption. One can imagine an experimental procedure, for instance using a mask, 
where graphene has regions exposed to molecular absorption separated from others 
where the concentration of adsorbates is small. This would produce a modulated electron 
concentration profile that would generate a modulated conductivity. 
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For illustrative purposes alone, we model the conductivity of graphene by 

cr(x) = ags{x) = ag[l — kcos{27[x/D)] , (1) 

where ag is the conductivity of homogeneous graphene and k is related to either the 
strain field field or the molecular doping concentration. Eq. ([T]) should be considered a 
toy model. 

The conductivity of graphene is a sum of two contributions: (i) a Drude term, 
describing intra-band processes and (ii) a term describing inter-band transitions. 
At zero temperature the optical conductivity has a simple analytical expression 
[22| [23t EU [251 ES]. The inter-band contribution has the form crj = cr'j + icr", with 

, / 1 hu — 2eF 1 hu + 2eF\ , . 

(T r = (Jo 1 H — arctan arctan , (21 

\ TT h'-y Tc nq J 

and 

1 {2ep + hujf + 
^^ = -"^2;^^" (2e^-M^ + ^V ^ ^ 

where cxo = vre^/ {2h). The Drude conductivity term is 

fTD = (Xo — , (4) 

vr Ai7 — inuj 

where 7 is the relaxation rate and ei;' > is the (local) Fermi level position with respect 
to the Dirac point. The total conductivity is 

ag = a' + ia" = a'j + laj + ao • (5) 




Figure 2. Optical conductivity of uniform of graphene: Drude (left) and inter-band 
(right) contributions. We assume ep = 0.45 eV and F — 2.6 mcV. The solid (dashed) 
line stands for the real (imaginary) part of the conductivity. 
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In Fig. |2] we plot the two contributions, Drude and inter-band, separately for a given 
value of e^;' and V = h'~^. For heavily doped graphene and for photon energies hujep < 0.5 
the optical response is dominated by the Drude term. Therefore, in what follows we 
assume 

c^g ~ o-D , (6) 

since we are interested in the regime of frequencies hw/ep ^ 0.5. For the frequency 
range of interest in this work (the THz spectral range) the above approximation gives 
accurate results. For the simulations given ahead we assume graphene sandwiched 
between two dielectrics of relative permittivities ei = 3 and €2 = 4, and an average 
Fermi level ep =0.45 meV, as in the experiments of Ref. \15l- In the same experiments 
the array of micro-ribbons has a period of 8 /im. In this work we choose D = 10 /xm. 



3. SPP dispersion 

Before solving the problem of ER scattering by a modulated conductivity, we overview 
some key concepts on surface plasmons-polaritons (SPPs) dispersion relation. SPPs are 
hybridized modes of electromagnetic radiation and free electrons of a conductor. They 
propagate along the interface of a dielectric and a conductor, and decay exponentially 
away from the interface. In the case of graphene, SPPs propagate along the graphene 
sheet. The SPPs have p— polarization (TM-wave) and their spectrum is given by 
1271 EH 122]: 

1 + — + = 0, (7) 

eiK2 uei 

and 

2 2 2/2 /o\ 

K^ = q -UJ /v^, (8) 

where m = 1,2 labels the media 1 (above) and 2 (below) the graphene sheet, 
Vm = \/^/{lJ'm^m) IS the velocity of light in medium m, q is the propagation wave 
number along the interface, and fim is the permeability of medium m (we shall assume 
that fim = f^o, the value in the vacuum). When the two media are equal, we obtain a 
simpler relation for the spectrum of the SPPs 

l + t^^^/^^^lJ^^=0. (9) 

The wave number gives the degree of localization of the SPP. This is defined as the 
ratio (assuming graphene in vacuum) [SD] 

— - 2af^ (10) 

CKi hhj 

where is the fine structure constant in free space and we have used Drude's formula 
for the conductivity, and took the limit F ^ hw. If we plug in typical numbers (see 
section ahead) we obtain 

uj 2 0.45 eV 



CKi 137 15 meV 



0.4, (11) 
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meaning a confinement size smaller than the wavelength of light in free space. Depending 
on the value of the Fermi energy and on the frequency of interest the degree of 
localization can be much smaller than 1. 

In general, Eqs. ([T]) and ^ cannot be solved analytically. However, approximating 
by its imaginary part (dispersive conductor approximation) 

ag ^ i^o^— , (12) 

and taking (for simplicity) ei = £2 = cq, the vacuum value for the permittivity, we 
obtain 



or 



f huj Y 



-ir) 



2 



1 + 



hcq Y 
afep) 



(13) 



(14) 



In the limit 2hcq/ {aftp) ^ 1 (non-retarded approximation) we obtain for the dispersion 
relation of the SPP the simple formula 

hjj = \j2afeFhcq . (15) 

In the limit g — i- we find that the SPP dispersion relation tends to w — )■ qc. For the 
case where ei 7^ €2, the non- retarded approximation yields 

hjj = \j2afeFhcq , (16) 

where 



^ A-nehc 

with e = (ei + e2)/2. We note that the behaviour of uj{q) as ~ cannot be trusted at 
small g, because for those values of q the non-retarded approximation is violated. In this 
case a full numerical solution of Eq. ([T]) is necessary. In what concerns the problem of 
ER scattering discussed ahead we will always be in the regime where the non-retarded 
approximation holds. 

In Figs. [3] and |4] we present the low momentum behaviour of the SPPs spectrum 
comparing the numerical solution of Eq. ([T]) for different levels of approximation to 
the optical conductivity and considering the role of the damping P = ^7. In Fig. |3] 
the effect of the damping is studied, keeping the full conductivity cxg = ctd- In this 
case the momentum q becomes a complex number, q = q' + iq", where q" describes the 
decay of the SPPs as it propagates in space along the graphene sheet. From this figure 
we see that the effect of the increase of the damping is two-fold: it shifts the SPPs 
dispersion relation toward higher energies and enhances the value of q" (as expected). If 
one wants to have long propagation lengths for the SPPs P must be as small as possible. 
In the same figure we also represent the light-hnes ficqj ^fe{ and hcq/ ^/ez which are the 
dispersion relations of a photon propagating in media 1 and 2, respectively. At grazing 



Light scattering by a medium with a spatially modulated optical conductivity: the case of grapheneJ 



ep.0.45 eV, e,.2, 6^-3 




0.01 0.02 0.03 0.04 0.05 0.06 



Figure 3. Plasmons-polaritons dispersion curves considering the effect of the 
broadening. The full conductivity ajj is included. The inset shows the dependence of 
the imaginary part of g", as function of the real part of g, q' . 

incidence it is possible to excite SPPs in graphene with a momentum q' corresponding to 
the point where the hght hne intercepts the dispersion curve of the SPPs. This occurs, 
unfortunately, at very low energies, v = u/{2n) < 1 THz, which corresponds to 4.2 
meV. 

r=2.6 meV; Ep.0.45 eV; e =1, e,=1 r=2.6 meV; Ep.0.45 eV; e =2. e,=3 




q' q' (1/nm) 

Figure 4. Plasmons-polaritons dispersion curves considering that the real part of the 
conductivity either finite or zero, and taking the cases of equal and different dielectrics 
sandwiching graphene. The straight lines marked ei and £2 are the light dispersion 
huj = hcq/y/el and Hu) = hcqj ^fl2 in media 1 and 2, respectively. 



In Fig. |4]we present the effect of neglecting the real part of the optical conductivity. 
The central feature is the vanishing of the dispersion curve for a finite value of g'. This 
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is a spurious result: when we use the full Drude conductivity the dispersion vanishes at 
zero q'. Changing the dielectric constants of the surrounding media changes the value 
of q' for the same frequency, as can be seen comparing the left and right panels of Fig. 

m 

When a periodic modulation of the conductivity exists, the SPPs develop a band 
structure depicted in Fig. |5| The black dashed lines represent the folding of the 



dispersion curve (16) into the first Brillouin zone, for vanishing modulation of the 
conductivity (k = 0). For finite modulation of a{x) {k ^ 0) gaps develop at the edges 
of the Brillouin zone. The band structure represented in Fig. [5] was computed in the 
non-retarded approximation, considering only the imaginary part of the conductivity 
(meaning that q = q' and q" = 0), and taking F = 0. As such, the behaviour of the first 
band close to g = is not accurate. In this approximation, there are no gaps at the 
center of the zone. (The approximations used allows to transform a non-linear eigen- 
problem into a linear one.) The region of the band structure between the two straight 
lines (see Fig. [5]) can be directly excited by light impinging on graphene without the 
aid of the ATR technique. If ER impinges at an angle 6, the momentum along the 
graphene sheet is fcisin^^ = w/fisin^^ and then a SPP with momentum q = cu/fisin^ 
can be created. This corresponds to a frequency of 



When the straight line (18) intercepts one of the upper bands in the Brillouin zone a 
SPP can be excited. If the modulation in the conductivity is weak (as, for example, is 
the case k = 0.3; see Fig. Isj) the energy of the excited SPP is roughly given by 



2afeFhc\ki s'm 6 — mG\ , (19) 

where G = 27i/D and m an integer. 

To conclude this section, we highlight the result we use later in the interpretation 



of ER scattering: for large q' the dispersion relation is approximately given by Eq. ( 16 ) 
and at high energies SPPs cannot be excited by shining light on homogeneous graphene, 
even at grazing incidence. When the optical conductivity is modulated periodically the 
system develops a band structure for the SPPs dispersion relation and their direct 
excitation becomes possible. 



4. Electromagnetic radiation scattering: formalism 

As discussed in Sec. [2] a conductive surface bearing a modulated conductivity supports 
SPPs which can be directly excited by light impinging on the conductor's surface, that 
is, without the need of a grating. We develop now the formalism needed to study, in a 
quantitative way, ER scattering by a modulated conductivity. 

Since the conductivity is a periodic function we can write it as a Fourier series 
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Figure 5. Band structure of the SPPs dispersion curves for tlie cosine modulation of 
the conductivity given by Eq. ([I]). Note that only the first six bands are shown. The 
light-cone, Eq. ( 18 ), at grazing incidence {6 = n/2) is represented by the dashed-dotted 
straight lines. 



(G = 27T/D): 



a[x) 



where 



AGmx^ 
e (TmG , 



m=—oo 



amG = (^{mG) = - I (T(x)e-*'"^^rfx . 



D 



(20) 



(21) 



Considering p— polarization, we write the electromagnetic fields as series of Bloch waves 



n=—oo 
oo 



izko sin9 ^ixko cos 6 j_ \ r^'^i^O sin 0—'t^G)z —q^^ \x\ 



56 , 

"T 2-^ '^z;k-nG'^ 
n=—oo 
oo 



A o izko sin 0ixko cose j_ \ ^ i{kosme-nG)z -qi^^\x\ /r)r)\ 

<Jl,\(ix;ili c -t- ^ ^x;k-nG^ ^ \'^'^) 



where A = 1, 2 labels the two cladding media. The relations between the amplitudes of 
the fields follow from Maxwell's equations: 

(i) incoming field {x < 0): 

sin 6 



cos 6 
ko cos 9 



(23) 
(24) 
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(ii) reflected field (x < 0): 



. ko sin 9 (i) 

^ (T) ^z;fcn ' 

./^OCl^ (1) 



(ill) transmitted field (x > 0): 

-,(2) 



_ /co sin 6^ (2) 



l(2) 



■ fM)^2(^ (2) 



The boundary condition i?^ = i?^ implies: 



,(1) 



,(2) 



The second boundary condition, — By = —fioaEl, implies 



t,(2) 



(25) 
(26) 

(27) 
(28) 

(29) 
(30) 

(31) 
(32) 



Using the relations between fields, Eqs. (24)- (28), the set of boundary conditions reduces 
to 



2ieiu 
kq cos U 

?^eS„ + ^eil+.5:a(^G)ea_ = A n^O. 



eic^ (2) 62^; (2) , ,^p^^^(J^e(2) 

(1) ^z-ko + (2) ^z;ko + ^ (^y^^)^z;k„ 
% % m 

eiw (2) e2Wj2) , „-\^^^_r-^j2) 

(1) ^z;kr, 



(1) ^z;kn ' (2) 

We note that q^^^ = —iko cos 6, with fco = (^/vi. We further define 

gi^) = ^(A;o sin^-nG)2-w2/^2^ 

for a positive argument of the square root (evanescent waves); if the argument of the 
square root in Eq. (35) is negative (propagating waves) qi!"^ is written as 



(33) 
(34) 

(35) 



g 



i\ u"^ /v\ — [kQsmO — nOy 



(36) 



The choice of sign for the square root is dictated by physical reasons: refiected waves 
for X < and transmitted waves for x > 0. The "pseudo-refiectance amplitude" of the 
mode n = is defined as 



rz;0 



^z;ko 



(2) 
^z;ko 



and that of a mode 7^ is given by 



(2) 



(37) 



(38) 
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The reflectance of tfie order n (for propagating modes) reads 

.(1) 2 



. kn COS 9 



(1) 



Tfie reflectance and tlie transmittance of tlie specular {n = 0) mode are given by 



and 



7^ 



r 



(2) 



£2 



COS 9 



ei Ves/ei - sin^ 9 



(39) 



(40) 



(41) 



respectively. The last expression is valid for e2 > ei or for 9 < 9c = arcsin(e2/ei), 
the critical angle for total reflection. In what follows only the the specular mode is 
propagating (all the remaining modes are evanescent). Then, the absorbance is simply 
given hjA=\ — lZ — T. 



5. Special limits 



Here we show that the linear system of Eqs. (33) and (34) reduces to well known 



formulae when the conductivity of graphene is either zero or finite and homogeneous. If 
the conductivity of graphene vanishes, we simply have an interface between two different 

= (n ^ 0)and 

-1 



(2) 

dielectrics. In this case Eqs. (|33p and (|34p give e\,.l^ 



(2) 



2 1 + 



cos 9 



ei \/l - tism^9/t2^ 
The reflectance amplitude follows from 

\/e2/ei cos 9 - A/T^^TTsii?^/^ 



1 - 



,(2) 
-'2;fco 



(42) 



(43) 



t a/ £2/ ei COS 6' + a/1 - ei sin^ 9/ €2 
which reproduces the well known result from elementary optics. 

Another particular limit is obtained when a{x) is flnite and homogeneous. In this 
case, only the Fourier component m = of a{x) exists. Thus, we obtain 



.(2) 



+ 



e2/ei 



e^-i I cos 6* A/e2/ei - sin^ 9 



+ 



eiVi 



cos 9 



(44) 



For 61 = €2 we obtain the well known result for the transmittance of graphene (for a 
TM wave) I26] 

n 2 



r 



2 + an cos^/eifi 



(45) 
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6. Electromagnetic radiation scattering: results 

When the optical conductivity of graphene has the form ([T]), the Fourier transform of 
s{x) reads: 

s{0) = 1 , (46) 



s{mG) 



(47) 



The linear system defined by Eqs. (33) and (34) is solved numerically and the sums 
over m are cutoff at m = — A^, . . . , 0, . . . , A^; the numerical solution rapidly converges 
with A^. 

Results for the specular reflectance, 71, and for the absorbance, A, are given in Fig. 
[6j As discussed in Sec. |3j it is not possible to excite SPPs in a homogeneous system 
because of the mismatch of the momentum of the incoming ER and that of SPP for 
the same frequency u. For a modulated conductivity the momentum of the SPPs is 



e =0.45 eV; r=0.6 meV; e =3; £,=4; k=0.3; D=10 [xm 
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Figure 6. Absorbance of a graphene sheet with a cosine-modulated conductivity. 
Reflectance (top) and absorbance (bottom) as function of the energy for different 
incoming angles. The parameters are ep = 0.452 eV, D = 10 ^m, t\ = 3, £2 = 4, 
/iF = 0.6 meV, and k — 0.3. The Brewster angle for the two dielectrics is 49.1°. 

conserved up to a reciprocal lattice vector mG, with m = ±1, ±2, . . ., that is 

ggpp = I /ci sin 6* — mG\ . (48) 

In this case, even for normal incidence, it is possible to excite SPPs. We stress that 
in the present case it is the periodicity of the conductivity and not an external grating 



leading to the validity of Eq. (48). The excitation of the SPPs at normal incidence is 
clearly seen in Fig. [6] In this figure both the reflectance and the absorbance are shown 
and a clear signature of SPPs excitation is present in both plots. The dashed black curve 
represents the behaviour of the system for a homogeneous conductivity and impinging 
ER at normal incidence; clearly the curve is featureless. For the inhomogeneous case a 
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large enhancement of the absorbance is seen around the energy given by Eq. (16) with 
q = 2ti/D. The position of the peak does not coincide exactly with the number given 
by Eq. ( |l6| ) because this equation is not sensitive to details of the band structure. From 
Fig. [5] we can see that the bands for m = ±1 at the zone center have an energy of about 
15 meV, coinciding with the energy for which the reflectance curve has a maximum for 
6' = 0. As the angle of the incident beam approaches the Brewster angle of the two 
media the reflectivity decreases substantially. The Brewster angle of the system is not 
given exactly by the well known formula Qb = arctan \ ft2ft\ due to the presence of 
graphene. For incidence angles above the Brewster angle the reflectance develops two 
dips and can be larger than it would be for 6* = (see curve for d = 80°). 

When the incoming beam deviates from normal incidence (i. e., 6 7^ 0), there is 
a peak splitting both in the reflectance and in the absorbance curves. We would like 
to understand in qualitative terms the behaviour of the peak slitting as function of the 
incoming angle 9. Using Eq. (|48|) in Eq. (|T6|) yields 



hcu = y2afeFhc\ki sin 6 — mG\ , (49) 

with ki = u/vi. Clearly, when m < and 9 increases the frequency shifts toward higher 
energies. On the other hand, when m > the energy decreases as 9 increases. This 
behavior can be understood from the analysis of Fig. [5j For ^ = the light line, Eq. 



(18), is vertical and touches the second band at the center of the zone (g = 0) where the 
branches associated with m = 1 and m = — 1 touch each other. As 9 grows, the slope 
of light line decreases and the branches with both m = 1 and m = — 1 are no longer 
degenerate. 



It is possible to obtain a simple analytical expression for uj(9) by solving Eq. (49) 
for hw. The final result is 

sm.9 1 



nuj=±— h — Vsin^6' + 4aci, (50) 



for m = —1 and m = 1, respectively, and the constants a and Ci are given by 

1 61 + 62 
" = ~A 



(51) 



In Fig. ^we plot the two branches of Eq. (50) and compare it with the position 



of the peaks for the absorbance obtained from Fig. |6l the agreement is qualitative. It 



cannot be quantitative because: (i) Eq. (50) is derived from a kinematic argument and 



therefore, misses the dependence on (ii) we have used the non-retarded approximation. 
However, for small k, the agreement is quite good. Indeed, if we had shifted the solid 
curves in Fig. [7]by a constant, they would overlap the points (solid squares) obtained 
from Fig. [g] In the bottom panel of Fig. [7]we plot \r±i\'^ for different ^ as a function 
of the energy. The energy of the peaks of \r±i\'^ coincides with that of peaks in the 
reflectance (and absorbance). This shows that the SPPs associated with the momenta 
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Figure 7. Position of the peak of the absorbance as function of the incoming angle. 
Top: dependence of the energy for maximum absorbance on the angle 9. The lines 
with squares are obtained from the bottom panel of Fig. |6] the solid lines are the 
two branches of Eq. ( 50 1 . The bottom panel shows the squared absolute value of the 
amplitudes of the modes associated with the SPPs of momentum ±G. The parameters 
are ep = 0.45 eV, D = 10 fim, ei = 3, e2 = 4, HT = 0.6 meV, and k — 0.3. 



-^2tt/D are the ones contributing to the anomahes in the reflectance and absorbance 
spectra. 

7. Conclusions 

We have shown that a modulated conductivity in an otherwise flat graphene sheet gives 
rise to a SPPs band structure with gaps opening at the edges of the Brillouin zone. 
Within the non-retarded approximation, no gaps open at the center of the zone. This 
band structure was later used to understand the qualitative behaviour of ER scattering 
from a flat graphene sheet with a modulated conductivity. The modulation allows for 
the SPPs be efficiently excited without using the ATR technique. The properties of the 
studied system resemble those of a photonic crystal. This analogy will be explored in a 
future work. 

If in addition to the modulated conductivity, corrugation is induced in the graphene 
sheet the reflectance and absorbance spectra are predicted to have a richer structure 
due to two competing mechanisms: the inhomogeneity and the presence of the groves. 
Finally we note that the case where graphene is cut into micro-ribbons is qualitatively 
different from the one discussed here, because the graphene sheet is, in that case, 
discontinuous. 



Light scattering by a medium with a spatially modulated optical conductivity: the case of graphenelh 
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